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Develop  the  MlNlMUF-3  model  to  provide  a simplified  prediction  of  hf  Maximum 
Usable  Frequency  (MUF)  suitable  for  use  on  small  mobile  propagation  forecast  (PROPHFT) 
terminals.  With  this  tool,  develop  a variety  of  new  forecast  applications  to  serve  the  hf 
surveillance  and  communications  community  where  the  use  of  large-scale  propagation 
programs  in  the  operational  environment  is  not  practical.  (The  tlexibility  of  these  new 
applications  to  serve  a variety  of  users  depends  to  a large  extent  on  the  confidence  of  the 
internal  prediction  algorithms.  Because  MlNIMUF-3  represents  such  a significant  degree  of 
simplification,  a primary  concern  is  the  accuracy  of  the  MUF  estimates.)  This  report 
presents;  (1)  the  MINIMUF  model;  (2)  results  of  verification  tests;  (3)  a specification  on 
its  capabilities  and  limitations;  and  (4)  applications. 


RESULTS 


The  MlNlMUF-3  hf  prediction  model  produces  a consistent  product  with  a nominal 
accuracy  between  3 and  5 MHz  rms  residual  error.  This  compares  favorably  with  the  abil- 
ities of  large-scale  hf  prediction  programs  which  require  sizable  computer  facilities. 


RECOMMENDATIONS 


Continue  work  to  further  refine  the  MINIMUF-3  model  to  include  E and  Sporadic  1. 
region  propagation.  Consider  small-scale  propagation  forecast  (PROPHET)  capabilities  for 
inclusion  into  new  hf  system  development.  Further  continue  new  application  development 
to  exploit  the  capabilities  of  simple  forecast  models  like  MINIMUF-3. 
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INTRODUCTION 


The  perforinance  of  any  system  operating  in  the  high-frequency  (hf)  spectrum  be- 
tween 2 and  32  MHz  is  subject  to  certain  physical  boundaries  on  the  propagation  of  its 
signal.  These  constraints  are  unique  and  definable  for  any  given  point  in  time  and  over  any 
path.  The  upper  boundary  is  called  the  maximum  usable  frequency  (MUF).  The  classical 
MUF  is  the  highest  frequency  that  can  propagate  by  a particular  mode  between  two  points 
by  ionospheric  refraction.  The  standard  MUF,  as  defined  by  the  CCIR  (ref  1 and  2),  is  an 
approximation  to  the  classical  MUF  obtained  by  combinine  data  from  vertical  ionospheric 
soundings  with  a simplified  analytical  solution  of  the  oblique  incidence  refraction  problem 
Wlien  observing  swept  frequency  oblique  incidence  sounding  data,  the  upper  frequency 
limit  is  referred  to  as  the  maximum  observed  frequency  (MOF).  The  MUF  (or  MOF)  is  a 
direct  function  of  ionospheric  electron  density  and  is  an  absolute  frequency  limit  in  that 
the  ionosphere  is  not  capable  of  supporting  higher  frequencies  for  that  path.  The  lower 
boundary  is  called  the  lowest  usable  (or  observed)  frequency  (LUF  or  LOF).  This  is  a 
function  of  absorption  occurring  during  daylight  hours  and  can  be  correlated  directly 
with  solar  zenith  angle.  The  band  of  frequencies  between  the  MUF  and  LUF  is  the 
propagation  bandwidth  or  the  spectrum  with  which  the  user  has  to  work.  The  cyclic 
nature  of  the  propagating  bandwidth  as  it  changes  from  night  to  day  is  shown  in  fig  1 
along  with  the  MOF  and  LOF  boundaries.  When  an  hf  system  operator  has  prior  knowl- 
edge of  the  changes  in  propagation  on  the  transmission  paths  he  wants  to  use,  he  can  adjust 
his  frequency  selection  to  remain  within  the  propagating  bandwidth.  When  he  has  no  such 
knowledge,  his  operation  is  vulnerable  to  these  propagation  changes. 

The  SOLRAD-Propagation  Forecast  (PROPHET)  concept  (ref  3)  is  based  on  pro- 
viding a near-real-time  assessment  of  the  hf  propagating  bandwidth  as  a function  of  time 
of  day,  season,  solar  activity,  and  transmission  path  geometry.  This  is  accomplished  by 
combining  established  and  new  hf  prediction  techniques  and  driving  the  resulting  forecast 
algorithms  with  geophysical  measurements  of  those  solar  emissions  which  have  an  impact 
on  propagation  parameters.  While  many  of  the  ideas  behind  PROPHET  have  been  tech- 
nically feasible  for  some  time,  recent  advances  in  mini-  and  microcomputer  technology 
have  made  the  use  of  these  concepts  also  economically  feasible.  The  initial  PROPHET 
concept  (described  in  ref  3)  uses  prestored  tables  generated  by  a standard  hf  prediction 
program  to  calculate  MUF  on  a large  computer.  Since  the  terminal  can  only  calculate 
paths  originating  or  terminating  at  the  location  for  which  the  tables  have  been  computed, 
the  MUF  algorithm  in  PROPHET  is  limited  to  the  fixed  paths  within  the  range  of  the 
tables  and  is  not  adaptable  for  mobile  operation. 


*■  International  Radio  Consultative  Committee,  13th  Plenary  Assembly,  Geneva  1974,  Recommendation 
373-3;  Definitions  of  Maximum  Transmission  Frequencies,  International  Telecommunications  Union, 
1975 

Bennington,  TW,  “Mow  Many  MUFs?,  Wireless  World,  v 65,  p 537-538,  December  1959 

Richter,  Dr  JH,  IJ  Rothmuller  and  RB  Rose,  PROPHET:  Real  Time  Propagation  Forecasting 
Terminal,  paper  presented  at  7th  Technical  Exchange  Conference,  El  Paso,  Texas,  30  November  1976. 
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Figure  1.  Cyclic  nature  of  propagating  bandwidth  as  it  changes  from  night  to  day. 


New  concepts  in  applying  PROPHFT  to  til' signal  security,  surveillance,  and  propa- 
gation tactics  required  a more  mobile  system.  To  accomplish  this,  the  MUF  calculation  had 
to  be  simplified  to  the  extent  that  it  could  be  run  on  new  microprocessor-based  systems. 


The  MINIMUF  model  was  developed  on  the  premise  that  the  majority  of  hf  predictions 
needed  for  the  employment  of  propagation  tactics  and  certain  hf  surveillance  missions  need 
only  be  first-order  approximations,  especially  in  a mobile  environment.  MlNlMUF-3  is  the 
current  stage  of  refinement  of  the  original  model. 


METHOD  OF  TESTING 


To  develop  a sufficient  knowledge  on  the  capabilities  of  MINIMUF-3,  two  verifica- 
tion approaches  were  considered.  The  first  was  to  compare  MlNlMUF-3  against  as  many 
swept  frequency  hf  oblique  sounder  data  as  could  be  assembled  in  a reasonable  time.  This, 
in  effect,  would  indicate  how  the  model  perfonns  when  correlated  with  “real-world”  pro- 
pagation. The  second  approach  was  to  compare  MINIMUF-3  results  against  large-scale  hf 
propagation  prediction  programs  such  as  ITS-78  (Red  Deck,  ref  4),  ITS-78  (Blue  Deck, 
ref  4),  and  ESSA-ITSA-1  (ref  5),  which  is  the  basis  for  the  NAVY  standard  NTP-6  hf 
predictions.  The  latter  type  of  comparison  would  provide  a relative  indication  of  the 
difference  between  a first-o’-der  approximation  and  a sophisticated  prediction  code.  While 
both  approaches  were  considered  useful,  heavy  emidiasis  was  placed  on  the  comparison 
with  oblique  sounder  data.  Since  the  large-scale  programs  were  considered  too  expensive 
to  rerun  for  the  project,  comparisons  with  the  results  of  other  prediction  programs  were 
performed  wliere  data  already  existed. 

OBLIQUE  SOUNDER  DATA  BASE  PREPARATION 

The  oblic|ue  sounder  data  base  assembled  for  the  MINIMUF-3  verification  was  de- 
rived from  a variety  of  sources  and  spans  the  period  between  19()0  and  1976.  Tliis  repre- 
sents more  than  one  complete  solar  sunspot  cycle  of  propagation  data.  Attempts  were 
made  to  make  the  data  base  as  diverse  as  possible  by  including  many  different  patli 
lengtiis.  orientations,  and  geographical  locutions. 

While  measurements  from  several  different  types  of  oblique  sounder  systems  were 
included,  the  majority  of  data  came  from  the  Navy  Tactical  Sounder  System  (NTSS).  To 
familiari/.e  the  reader  with  oblique  sounding  principles,  a brief  description  of  the  NTSS 
follows. 


Idiviromncntal  Sciences  .Services  Administration  technical  report  KRL-I IO-ITS-78, 
Predicting  Long  Tenn  Operational  Parameters  of  High  Frequency  Sky  Wave 
Telecommunication  Systems,  by  AF  Barghausen,  JW  Finney,  LL  Proctor  and 
LD  Schultz,  May  1969. 


Fnvironmcntal  Sciences  Services  Administration  technical  report  IFR-l-lTSA-l 
Predicting  Statistical  Performance  Indexes  for  High  Frequency  Ionospheric 
Telecommunication  Systems,  by  DL  Lucas,  ('.W  Haydon,  et  al,  August  1966 
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NAVY  [ ACTICAL  SOUNDhR  SYSTHM  (NTSS) 


The  Navy’s  oblique  ionospheric  sounder  system  consists  of  several  shore-based  | 

sounder  transmitters  and  a miber  of  sounder  receivers.  AN/FPT-1 1 (XN-1 ) sounder  | 

transmitters  are  installed  at  selected  Naval  Communications  Stations,  research  installations,  t 

and  aboard  shi()s.  . 

AN/TPT-1  1 Transmitter.  Once  each  minute,  the  FPT-1  1 transmitter  sequentially 
transmits  a pulse  on  each  of  80  discrete  channels  between  2 and  32  MH/;  the  total  scan 
consists  ot  160  pulses  and  has  a scan  time  of  16  seconds.  The  frequency  range  is  divided 
into  tour  active  bands,  with  20  channels  linearly  spaced  within  each  band;  the  spacing 
between  adjacent  channels  is  doubled,  proceeding  from  one  band  to  the  next  band.  The 
spacings  arc  100  kll/  in  the  2-to^-MHz  range  (Band  A);  200  kHz  in  the  4-to-8-MHz  range 
(Band  B);  400  kHz  in  the  8-to-16-MHz  range  (Band  C);  and  800  kHz  in  the  16-to-32-MHz, 
range  (Band  D).  These  spacings  provide  a 57e  frequency  resolution  at  the  low  frequency 
limit  ot  each  band  and  a 2.59f  frequency  resolution  at  the  upper  frequency  limit  of  each 
band.  Tabic  I shows  the  exact  channels  used  by  the  FPT-1 1.  Each  frequency  transmitted 
consists  ot  two  pulses  separated  by  50  milliseconds.  Each  2.6  ms  pulse  is  composed  of  a 
series  ot  13  subinilses,  biphase-coded  in  a Barker-code  sequence.  The  use  of  Barker 

coded  pulse  compression  increases  the  effective  transmitter  power,  retains  necessary  } 

pulse  resolution,  and  increases  the  received  signal-to-noise  ratio  by  13;  I . To  eliminate 
interterence  between  transmitters,  starting  times  of  the  1 6-second  scan  are  staggered, 
each  transmitter  site  being  assigned  a precise  time.  The  transmitter  power  is  normally 
30  kW  (PF.P).  The  antenna  must  be  wideband,  50-  to  72-ohm  impedance,  and  have  a 
VSWR  of  less  than  3: 1 over  the  band  of  2 to  32  MHz. 

; AN  UPR-2  Receiver.  The  AN/UPR-2  receiver  sequentially  processes  the  pulse-train 

i input  by  starting  the  gated  receiver  scan  at  the  time  of  the  transmission.  This  is  accomp- 

i lished  by  synchronizing  to  a common  timing  source  (ie,  WWV)  and  maintaining  an  accurate 

r time-base  generator  in  the  receiver.  Signal  processing  is  required  since  each  sounder  signal 

I is  composed  ot  a series  ot  13  subpulses;  each  subpulse  compression  is  such  that  the  signal 

pulse  voltages  add  arithmetically  but  the  noise  does  not.  The  improvement  in  signal-to- 
I noise  ratio  (process  gain ) is  the  voltage  ratio  of  N : 1 , where  N is  the  number  of  RF  pulses 

in  the  seciuence.  For  the  13-bit  code,  the  process  gain  is  1 1 dB.  (The  Barker  code  is  used 
\ because,  after  processing,  it  produces  a video  pulse  having  a 13: 1 peak-to-sidelobe  ratio, 

I which  was  found  to  be  optimum  tor  13-bit  codes.)  The  resultant  or  summed  pulses  are 

I then  fed  to  a storage  oscilloscope,  which  is  used  to  present  this  pulse  information  in  video 

! form.  A storage  oscilloscope  incor|iorated  in  the  sounder  receiver  displays  and  retains 

; each  scan. 
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Table  1.  AN/F'PT-1 1 Step  Sounder  Frequency  Channels  (2-32  MHz). 


BAND  A 

(MHz) 

BANDB 

(MHz) 

BAND  C 

(MHz) 

BAND  D 

(MHz) 

2.075 

4.14 

8.3 

16.6 

2.175 

4.35 

8.7 

17.4 

2.275 

4.55 

9.1 

18.2 

2.375 

4.75 

9.5 

19.0 

2.475 

4.95 

9.9 

19.8 

2.575 

5.15 

10.3 

20.6 

2.675 

5.35 

10.7 

21.4 

2.775 

5.55 

1 1.1 

22.2 

2.875 

5.75 

11.5 

23.0 

2.975 

5.95 

1 1.9 

23.8 

3.075 

6.15 

12.3 

24.6 

3.175 

6.35 

12.7 

25.4 

3.275 

6.55 

13.1 

26.2 

3.375 

6.75 

13.5 

27.0 

3.475 

6.95 

13.9 

27.8 

3.575 

7.15 

14.3 

28.6 

3.675 

7.35 

14.7 

29.4 

3.775 

7.55 

15.1 

30.2 

3.875 

7.75 

15.5 

31.0 

3.975 

7.95 

15.9 

31.8 

i 

This  display,  shown  in  fig  2,  indicates  which  frequencies  are  propagating  [lowest  observed 
frequency  (LOF)  to  maximum  observed  frequency  (MOF)]  and  which  frequencies  will  be 
highly  distorted  because  of  multipath  time  delay.  An  additional  data  output  from  the 
standard  UPR-2  is  a strip  chart  record,  which  is  produced  concurrently  with  each  ionogram 
record.  These  outputs  provide  a permanent  record  of  the  daily  variations  of  the  scanned 
spectrum  between  2 and  32  MHz. 
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l igurc  2.  loiwgram  of  oblique  sountlcr  transmission  from  Naval  Communications 
Station,  Hawaii. 

IK  Digital  Recorder  (HFDR).  To  .supplement  this  capability,  NOSC  developed  a 
method  ot  digitizing  the  video  output  signal  and  recording  it  on  magnetic  tape.  This 
inethod  evolved  into  the  design,  construction,  and  implementation  of  the  hf  digital  recorder 
(HFDR)  shown  in  tig  3. 
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The  hf  sounder  digital  data  recording  system  is  a self-contained  unit  consisting  of 
two  primary  subsystems: 

a.  The  analog  to  digital  (A/D)  and  timing  interface  (DIGITIZER);  and 

b.  The  magnetic  tape  transport. 

With  minor  modification  to  the  AN/UPR-2,  the  HFDR  operates  concurrently  with 
the  AN/UPR-2  receiver  and  in  no  way  affects  normal  operation.  Hence,  with  the  HFDR-II- 
equipped  sounding  receiver,  all  amplitude,  time  delay,  and  frequency  information  is  re- 
corded once  every  minute,  24  hours  a day. 

OTHER  SOUNDER  SYSTEMS 

Data  collected  prior  to  1968  were  measured  on  a variety  of  systems.  These  included 
Graiiger  900  systems.  Modified  Granger  1900  systems.  Modified  C3  ionosondes,  and  special- 
ly developed  systems.  In  general,  the  frequency  steps  used  by  the  Granger  systems  are 
compatible  with  the  NTSS  receivers  and  in  some  instances  data  were  acquired  by  means  of 
a Granger  transmitter  and  UPR-2  receiver. 

DATA  CATEGORIZATION 

Because  of  the  variety  of  measurement  systems  and  recording  methods,  some  com- 
mon MUF  parameter  had  to  be  established  for  the  comparison  with  MlNlMUF-3.  Since 
most  of  the  data  could  be  reduced  to  an  hourly  mean  MUF  as  a function  of  month  and 
year  (hence  sunspot  number),  the  path-mouth  was  established  as  the  basic  data  block  for 
comparison  purposes.  Each  path-month  consists  of  24  hourly  mean  MUFs  for  a given 
transmission  path. 

Tlte  source  of  the  oblique  sounder  data  was  also  deemed  important  because  it  in- 
fiuences  the  statistical  significance  of  a given  path-month.  The  overall  sounder  data  were 
categorized  into  three  sources: 

a.  NTSS-HFDR 

b.  NTSS-Strip  Chart 

c.  Non-NTSS. 

A path-month  MUF  curve  from  the  NTSS-HFDR  system  is  generally  the  product 
of  approximately  40000  digitally  processed  measurements.  The  resolution  of  the  NTSS- 
strip  chart  system  limits  this  to  about  2880  hand-scaled  data  points  per  path-month.  Be- 
cause the  documentation  on  the  non-NTSS  Systems  was  insufficient  with  respect  to  data 
reduction  techniciues,  a path-month  is  assumed  to  consist  of  30  points  per  hour  or  720 
data  points.  For  these  reasons,  the  investigators  in  this  study  placed  the  greatest  con- 
fidence in  the  NTSS-HFDR  data. 


DESCRIPTION  OF  OVERALL  SOUNDER  DATA  BASE 

The  final  obliiiue  sounder  data  set  consists  of  196  path-months  (or  16.33  path-years) 
of  mean  hourly  MUF  values  derived  from  23  different  hf  transmission  paths.  The  longest 
path  is  7808  km  and  the  shortest  path  is  802  km.  The  set  contains  a cross  section  of 
transmission  paths  including  midlatitude,  transauroral,  transequatorial.  all  seasons,  and  all 
solar  sunspot  numbers  (SSN ).  No  single  path  represents  more  than  1 2%  of  the  total  data 
set.  Table  2 summarizes  the  basis  against  which  the  MINIMUE'  model  was  compared. 


Table  2.  ill  Propagation  Obliiiuc  Sounder  Data  Base. 


(14)  Honolulu,  Hawaii  21.42”N  158.15  W 4200  - 22  1 1 1968-70  87-103  NOSC  file  NTSS  (ST) 

to  Corona.  CA  33.86°N  117.55°W 


Table  2.  (Continued) 
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While  sounder  MUF  data  were  available  in  the  NOSC  data  file,  it  was  believed 
that  the  paths  listed  in  table  2 represented  a sufficient  cross  section  of  propagation  condi- 
tions to  test  MlNlMUF-3  thoroughly. 

Because  the  data  came  from  a variety  of  sources,  a standard  format  was  established 
and  a master  data  tape  was  constructed  listing  the  hourly  (UT)  mean  MUFs  and  all  pertinent 
path  parameters.*  Once  constructed,  the  desired  comparisons  to  the  MINlMUF-3  model 
were  a relatively  straightforward  process  using  the  NOSC  SEL810B  computer  system. 


* A copy  of  this  tape  can  be  made  available  upon  request  through  Commander,  NOSC. 
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MINIMUF  MODEL  DEVELOPMENT 
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[ 

I 

( 
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THEORETICAL  BACKGROUND 

A central  task  olTong-term  statistical  hf  propagation  forecasting  is  the  prediction  of 
maximum  usable  freciuency  (MUF-  ).  The  MUF,  in  turn,  is  principally  controlled  by  the 
critical  frequency  of  the  I-’-i  layer  of  the  ionosphere,  10^2’  success  in  predicting 

this  quantity  that  primarily  determines  the  accuracy  of  the  MUF  forecast.  Unlike  the  F and 
F]  layers,  which  can  be  modeled  quite  well  as  a function  of  a single  parameter- cos  x (the 
cosine  of  the  instantaneous  solar  zenith  angle)  proportional  to  the  solar  intensity  the  physics  j 

of  the  Ft  layer  is  generally  believed  to  involve  an  interaction  of  photochemical  and  trans- 
port processes  sufficiently  complex  that  diurnal,  seasonal,  and  geographical  I'^Ft  variations 
cannot  be  simply  accommodated  through  the  corresponding  variations  in  cos  x-  Indeed, 
one  even  speaks  of  Ft  layer  “anomalies”  when  comparing  observed  fo^2  '''***’  expectations 
based  on  the  instantaneous  cos  x (ref  1 2 and  13).  For  example,  f^FT  can  be  higher  at 
midday  in  winter  than  in  summer  (“seasonal”  anomaly),  and  on  a given  day  can  peak  in 
late  afternoon  rather  than  at  midday  (“diurnal”  anomaly).  ; 

Therefore,  while  1^1-2  cannot  be  modeled  as  a function  of  the  instantaneous 
cos  X,  the  possibility  remains  that  it  could  be  modeled  as  the  response  of  a dynamic  system  } 

“driven”  by  a function  of  cos  x.  lixamination  of  the  shape  of  observed  f^jFT  diurnal  pro- 
files, for  example,  suggests  that  a simple  relaxation  model,  according  to  which  I'^Ft 
represents  a lagged  response  to  the  instantaneous  solar  intensity,  may  be  useful  as  a first 
approximation.  Allowing  the  lag  time  constant  to  be  long  (~  10  hours)  in  summer  and 

short  (~  1 hour)  in  winter  at  middle  and  equatorial  latitudes  could  then  at  least  partially  ! 

reproduce  both  the  seasonal  and  diurnal  anomalies.  ■ 

Accordingly,  a semiempirical  model  for  fQF2  has  been  developed  based  on  the  ‘ 

analogy  to  a single-lag  linear  system  (eg,  an  RC  circuit)  driven  by  a forcing  function  j 

proportional  to  the  instantaneous  cos  x.  Further  simplifying  assumptions  of  the  model  j 

are  as  follows:  1 

1 . The  lag  time  constant  during  the  day  is  a simple  monotonic  function  of  the 

midday  solar  zenith  angle.  :| 

2.  The  time  constant  at  night  is  a constant  (2  hours)  independent  of  season 

or  geographical  location.  j 

As  with  other  semiempirical  models  of  complex  geophysical  processes,  no  attempt  is  made  S 

to  justify  the  model  in  terms  of  the  underlying  physical  mechanisms.  Rather,  the  model  | 

serves  to  provide  a mathematical  framework  for  force-fitting  to  empirical  data.  Of  course.  j 

if  the  model  is  successful  in  fitting  a large  data  base  with  reasonable  accuracy  and  rela-  \ 

lively  few  adjustable  constants,  the  physical  reality  of  the  assumed  relaxation  process  J 

gains  credibility  and  may  guide  the  understanding  of  the  underlying  mechanisms.  ^ 


12  JA  Ratcliffc  and  K Wcekcs,  in  Physics  of  the  Upper  Atmosphere,  J A Ratcliffe  (cd ),  p 427  432,  | 

Academic  Press,  New  York,  I960. 

13  Rishbeth,  II,  A review  of  Ionospheric  I-  Region  Theory,  Proc  IliFH,  v .S.S,  p 16  3.S,  January  1967. 
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A key  toaturo  ot  the  proposed  IqI'i  model  is  that  seasonal  and  geographical  vari- 
ations ot  the  predicted  lyFi  arise  only  from  the  corresponding  variations  in  the  midday 
solar  zenith  angle,  in  marked  contrast  to  the  customary  procedure  of  numerically  mapping 
*0^  2 tilting  appropriate  mathematical  functions  to  observed  ionospheric  sounding  data 
collected  from  a worldwide  net  of  vertical  sounders  (ref  14,  15,  and  16).  Furthermore,  by 
making  simple  analytical  approximations  to  the  dynamic  solutions  of  the  model  (ie,  the 
diurnal  response  (unction),  we  end  up  with  a simple  closed-form  expression  for  foF2  as  a 
tunction  ot  midday  solar  zenith  angle,  sunspot  number,  and  time  relative  to  local  sunrise 
and  sunset.  By  appending  simple  approximations  for  the  M factor  (ie,  MUF/f^F-i)  and 
lor  solar  zenith  angle  as  a function  of  location  and  time,  we  arrive  at  a model  for  MUF 
which  is  suttlciently  compact  to  be  coded  for  computation  on  a minicomputer  or  desk- 
top programmable  scientific  calculator. 


MATHEMATICAL  MODELS 


Let  cos  X be  the  cosine  of  the  actual  instantaneous  solar  zenith  angle  at  an  iono- 
spheric height  representative  of  the  F->  layer,  and  let  cos  X<.ff  the  “effective”  value  of 
this  quantity  taking  relaxation  processes  into  account.  The  basic  assumed  form  for  f^F-i  is 


'of':  - ( ‘■'■f^/Ro)  (V'^osxeff  ) 


(1) 


where  R - sunspot  number  and  Rq,  Aq,  and  A|  are  constants  independent  of  geographic 
location  and  time*.  To  model  cos  x^ff'  we  begin  by  constructing  the  quantity: 


^^'osXc-ff 


SUNSET 


•exp 


(2) 


!4  International  Radio  Consultative  Committee,  Oslo,  1966,  Report  340;CCIR  Atlas  of  Ionospheric 
Characteristics,  International  Telecommunications  Union,  Geneva,  1975  . 

15  International  Radio  Consultative  Committee,  New  Delhi,  1974,  Supplement  No  1 to  Report  340; 

( ( IR  Atlas  o(  Ionospheric  Characteristics,  International  Telecommunications  Union,  Geneva,  1971. 

16  International  Radio  Consultative  Committee,  Geneva,  1974,  Supplement  No  2 to  Report  340; 

( ( IR  Atlas  ol  Ionospheric  Characteristics,  International  Telecommunications  Union,  Geneva,  1975. 

* This  lorm  was  motivated  by  the  following  physical  considerations:  ( 1 ) at  the  peak  of  a Chapman 

layer,  the  solar-induced  electron  production  rate  is  acos  x;  (2)  assuming  the  net  effect  of  the  electron- 
loss  chemistry  to  be  recombinationlike,  the  solar-induced  increase  in  electron  density  over  ambient 
is  cv  (production  rate)^;  (3)  for  a collisionless  plasma,  critical  frequency  is  afeleciron  density)^. 
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In  equation  2,  rjq  is  a nighttime  relaxation  time  which  is  taken  to  be  a constant  independent 
ot  season  and  geographical  location.  T - is  the  elapsed  time  in  hours  since  sunset. 

During  the  day,  a different  relaxation  time,  tq,  is  assumed.  In  contrast  to  rj^j, 

markedly  dependent  on  location  and  position,  albeit  in  an  analytically  simple  fashion. 
Specitically,  we  assume  the  daytime  relaxation  time  to  be  a function  of  the  actual  noon- 
time solar  zenith  angle  (cos  XjsjoON 


= ^o  ( 


xnoonJ 


where  and  Pt  are  again  constants  independent  of  season  or  location.  Note  that  during 
summer  at  equatorial  and  moderate  latitudes,  Tp  -*•  Tq  , whereas  in  the  winter  . 

Tlie  time  dependence  of  cos  during  the  day  is  approximated  as  follows,  f'irst, 
we  simplily  the  time  dependence  of  the  actual  cos  x in  tenns  of  its  value  at  noon, 
cos  XnoON  ’ 


cos  X ==  cos  XNOON  ' 


^ ^'^■t'dawn^ 


where  AT  is  the  daytime  duration  given  by 


- "^SUNSET -T'daWN 


and  noon  occurs  when  T - + AT/2.  Next,  cos  x^ff  is  assumed  to  represent  the 

response  of  a linear  first-order  system  “driven”  by  the  actual  cos 


’’D  Hf  ‘‘^°*>:eff>  cosXeff=cosx. 


Using  equation  4 for  the  right  hand  side  of  equation  6 and  integrating,  we  obtain  for  the 
daytime  cos  x^;}-f: 


\ _ cos  XNOON  ( / '"^-TdaWN' 


where 
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(9) 


uiul 


a 


AT 


I’valiiating  cciuation  7 at  sunset  = AT),  we  find 

/ \ xnoon 


jsUNShT 


+ P 


I 


( 

r / \ 1 

, /-at\ 

1 

1 + exp  1 ) 

/ 

t 

\ /J 

) 

(10) 


wliieli  eonipletes  the  specitleation  of  (qI  t nt  night  via  equations  1 and  2. 

F inally,  to  avoid  unphysieal  diseontinuities  in  cos  x,.ff  just  after  suniise,  we  do  not 
allow  cos  ^ to  fall  below  its  value  just  prior  to  sunrise  as  computed  from  equation 

2.  Thus  we  impose  the  condition 


Ihe  basic  f^jIS  model  is  therefore  given  in  equations  1, 2,  and  10  for  nighttime 
(ie,  I < T<T[)awn  * I,  7,  and  II  for  daytime 

‘Tdawn^^'^t^sunsf.t*- 


i;giiATioN  OF-  TiMi; 

In  this  section,  we  present  simple  analytical  approximations  for  the  times  of  local 
iU)on,  sunrise,  and  sunset,  and  for  the  noon  value  of  the  solar  zenith  angle.  To  an  accept- 
able degree  of  aceuracy,  we  find  for  the  Universal  Time  of  local  noon  the  approximation; 


'noon  =7^  +12+0.12 


sin  Y I + 1 .2  sin  (2Y  j ) 


mod  (24) 


where 

W = longitude  west  ofCireenwich  [0,360°  | 
Y,  = 0.0172  (I)+10) 

I)  = 30.4  (M 1-1 ) + 1)1 
1)1  = day  [1.311 

Ml  = month  [ I,I2[. 


(12) 
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In  terms  of  the  subsidiary  variable  defined  by 


= 0.409  cos  Y ] [radians] 

and  L = north  latitude  in  degrees  [-90,  90] 
we  have 


i 

a \ 

1 

cos  — 

C“7 

1 + yo 

2 \ 

f - 

cos  XnoON 

The  duration  of  the  daytime,  AT,  is  approximated  by 

/ -0.26  + sin  yT  sin  L \ 


24 

AT  = — arc  cos 

TT 


cos  yo  cos  L 


[hours] 


(13) 


(14) 


(15) 


where  the  factor -0.26  approximately  represents  the  difference  between  sunrise  (or  sunset) 
at  the  surface  of  the  earth  and  at  Ft  layer  heights.  From  equation  15  there  follows 

T'dAWN  = ’rNOON"^T/2  mod  (24)  (16) 

^SUNSET  " ^NOON 
where  TjqQQjq  is  given  in  equation  1 2. 


CONTROL  POINTS 


The  equations  of  the  foregoing  sections  yield  foF2  ^ specified  (north)  latitude,  L, 
and  (west)  longitude.  W.  In  an  actual  application,  the  latitude  and  longitude  of  the 
receiver  and  transmitter  are  given  quantities  and  L and  W are  to  be  evaluated  at  specific 
“control  points’’  along  the  great  circle  propagation  path. 

From  spherical  trigonometry,  one  obtains  the  following  expressions  for  the  great 
circle  arc  length,  i//  (in  radians),  connecting  two  points  defined  by  the  (latitude,  longitude) 
pairs  (LpW | ) and  (Lt.Wt): 


= arc  cos 


sin  L|  sin  Lt  + cos  Lj  cos  Lt  cos  (Wi-W] ) 


(17) 


If  one  travels  a fraction,  K,  of  the  distance  from  point  1 to  point  2,  the  north  latitude,  L|^. 
and  west  longitude.  Wj,^,  of  this  location  are  detemiined  from 


arc  sin  j 

[sin  i^(l-K)  sin  L|  + sin  ((i'K)  sin  Lt  | 

sin  i//  1 

(18) 
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and 


Wj^  = arc  cos 


I cos  L[  cos  W I sin  i//(l-K)  + cos  Lt  cos  Wt  sin  (i//K ) ) 

^ = J 

cos  L|^  sin  \jj  I 


09) 


rile  control  point  locations  are  detemiined  in  accordance  with  conventional  MUF 
estimation  procedures.  Thus  it'  the  range,  D,  from  transmitter  to  receiver  is  less  than 
4000  km.  the  control  point  is  at  midpath  (K=V2).  If  D > 4000  km,  two  control  points 
are  considered,  respectively  located  2000  km  in  from  either  tenninus.*  In  this  latter  case, 
the  procedure  we  adopt  is  to  compute  a MUF  for  each  of  two  fictitious  single-hop  paths  - 
ot  length  D/2,  whose  midpaths  lie  at  the  control  points  just  defined.  The  MUF  for  the 
actual  total  path  is  then  taken  as  the  lesser  of  the  two  single-hop  MUFs.  Thus  for  any  D, 
the  problem  reduces  to  the  estimation  of  MUFs  for  a single-hop  path. 


M FACTOR 

Again  following  conventional  procedure,  the  MUF  for  any  single-hop  path  is 
written  as  the  product: 


MUF^U'oFt) 


mid-path 


•M 


<20) 


where  the  so-called  “M  factor”  is  a function  of  pathlength  and  layer  structure -particularly 

the  height  of  the  layer  maximum.  In  the  present  model,  (f^F-))  is  computed  at  a 

^ mid-path 

control  point  as  fully  described  above.  We  therefore  burden  M with  accommodating  the 
important  “real-world”  complexities  of  hf  propagation  which  have  been  glossed  over  in 
the  extremely  simple  models  we  have  adopted  for  f^F^. 

Consider  a single-hop  path  of  great-circle  arc  length  \p  (radians)  whose  midpoint 
lies  at  latitude  L^.  and  whose  termini  lie  at  latitudes  L]  and  L-(.  Let  the  duration  of 
daylight  at  midpath  be  AT  hours.  We  approximate  the  M factor  in  the  following  form 


M 


= |l+2.5[sin(2.5<//)l^/-|  •G,(AT)-G2(Lo)-G3(L,,L2). 


(21) 


7 he  first  factor  of  equation  2 1 embodies  the  range-dependence  of  the  M factor  and 
was  obtained  by  curve-fitting  to  exact  results  for  a parabolic  layer  of  290-km  height,  with 
a ratio  ot  seniithickness/base  height  = 0.4.  The  G|  factor  gives  recognition  to  the 
increase  in  ['2  layer  heights  observed  at  high  (northern)  latitudes  during  the  summer; 
ie,  at  or  near  “midnight  sun”  conditions.  We  specifically  choose 


G]  = I -O.lexp  ((AT-24)/31  (22) 

which  has  the  effect  of  a lO'/f  reduction  in  MUF  under  full  midnight  sun  conditions,  with 
(i|  recovering  rapidly  'vl  as  the  (midpath)  latitude  moves  toward  the  equator. 

* I km  = 0.6  statute  mile  or  0.5  nautical  mile. 
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C'l-)  is  designed  to  produce  further,  seasonally  independent,  reductions  in  MUf'  for 
high-latitude  oaths.  Since  propagation  data  point  to  lairly  abrupt  onset  ol  this  reduction 
for  I latitudes  | > 4.*)°,  we  choose  the  step  fomi : 

1 if  1 latitude  l<  45°  ^ 

" 0.8  if  1 latitude  45° 

ie,  a 20%  MUF  reduetion. 

Finally,  (jj  is  a correction  factor  for  transequatorial  paths  to  take  approximate 
account  of  the  well-known  MUF  increases  on  such  paths.  Applying  another  209?  correction 
in  step-wise  fonn: 


1 if  Li  , L-)  of  same  sign 

Gi  = 

1 .2  if  L I . L-i  ot  opposite  sign 


FIT  TO  OBLIQUE  SOUNDER  DATA 

The  as  yet  unspecified  constants  R^,  Aq,  A|  (eq  1 );  T]sj  (eq  2);  and  r^.  Pt  (eq  3 ); 
together  with  the  correction  factors  (i] . Gt  and  G3,  were  determined  by  iteratively  ad- 
justing the  model  to  portions  of  the  data  base  described  earlier.  The  specific  procedure 
followed  was  to  first  identify  a limited  subset  spanning  the  range  of  path  types  (ie,  long 
and  short  paths,  high  and  moderate  latitudes,  north-south  and  east-west,  etc).  Then  tor 
this  set,  representative  path  months  were  chosen  and  the  observed  diurnal  MUF  variations 
were  simultaneously  displayed  on  a high-resolution  CRT,  using  a minicomputer-based 
Megatek  graphics  display  system.  The  model  was  implemented  in  BASIC  in  the  mini- 
computer so  that  direct  simultaneous  visual  comparisons  could  be  made  for  the  whole 
subset  (36  path  months)  between  the  observed  and  predicted  diurnal  MUF  variations. 

Since  computation  and  display  of  the  864  MUFs  (ie,  36  path  months  X 24  hours)  took 
about  10  minutes,  interactive  “tuning”  of  the  model  was  quite  feasible  and  efficient, 

Tlie  final  parameter  values  obtained  in  this  way  are 

Ro  = 250  (SSN) 

Ay  = 6 (MH/'^'-) 

A]  = 58  (MH?/^) 

T,q  = 2 (h) 

Tq  = 9.7  (h) 

Pi  = 9.6  (exponent-no  dimensions). 

It  is  evident  from  the  “roundness”  of  these  numbers  that  further  fine-tuning  of 
the  model  on  an  extended  data  base  is  still  possible.  For  the  present  verification  purposes, 
the  model  was  frozen  at  these  parameter  values  and  was  designated  MINIMUF'-3.  BASIC 
and  FORTRAN  versions  are  given  in  the  appendices. 


MlNIMUF-3  VERIFICATION  RESULTS 


riiis  suction  will  present  tlie  results  of  the  MINIMUF'3  verification  tests.  The 
objective  is  to  provide  the  reader  with  a clear  understanding  of  the  capabilities  and  limi- 
tations of  the  new  MUF  prediction  model. 

COMPARISON  CRITERIA 

The  primary  comparison  parameter  used  throughout  the  verification  test  was  RMS 
error  (F,.,„j^»  between  observed  oblii|ue  sounder  MOF's  and  MUFs  calculated  by  MINlMUF-3. 
This  was  calculated  by 

* rms  ~ I 

where  m^^  = observed  hourly  mean  MUF'  for  a given  path  month 

uij^.  = calculated  MUF'  for  that  hour,  path,  midmonth  day,  and  ^SN 
n = number  of  data  compared. 

file  nils  error  was  used  because  it  would  show  the  ability  of  MlNlMUF-3  to  track  MUF 
variations  along  with  providing  error  bounds.  In  addition,  the  use  of  nns  error  would 
allow  direct  coni|ntrison  of  results  to  previous  studies  where  large-scale  hf  prediction 
programs  were  compared  to  some  of  the  same  hf  obliciue  sounder  data  used. 

The  iirincipal  method  used  to  assess  the  accuracy  of  MINIMUF-3  was  to  derive 
the  nns  error  for  each  jnitli  month  in  the  overall  data  base  and  assem^'le  these  in  a variety 
of  categories  to  determine  the  strengths  and  weakness  of  the  model.  These  comparisons 
are  |iresented  in  the  following  section. 


RESULTS  OF  MINIMUF-3  VERSUS  SOUNDER  DATA  COMPARISON 

I.  ALL.  ( ASF. S 

Hie  obliciue  sounder  data  base  consisted  of  196  path  months  (16.3  path  years)  of 
observed  MUFs  taken  over  23  different  hf  transmission  ( ■'tlis.  The  shortest  path  was  802 
km  and  the  longest  was  7808  km.*  Table  3 lists  each  of  the  paths  and  the  total  mis  error 
tor  that  path. 

For  the  entire  Uhi  path  months  (4704  hourly  points),  the  overall  rtiis  error  of 
.MINIMl)F-3  was  3.83  Mil/..  Within  this  error,  the  accuracy  of  61%  of  the  sample  was 
eiiual  to  better  than  the  design  goal  of  4.0  MHz.  The  remaining  39%  of  the  tests  produced 
mis  errors  between  4.05  and  5.82  MHz. 


* 1 kill  - O.Ci  sluline  mile  or  O..*!  iiaiitieal  mile. 
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Table  3.  MINlMUF-3  rnis  Error  by  Sounder  Path 
(Total  Sample  = 196  Path  months  of  Data) 


No. 

Transmission  Path 

Type 

Path 

Length, 

km  ** 

Solar 

Sunspot 

No 

9c  of 
Sample 

MINIMUF-3 
rms  Error, 
MHz 

1 

Palo  Alto,  CA,  to  Fairbanks,  Alaska 

H,N/S 

3503 

17 

1 

2.03 

2 

Andoya,  Norway,  to  Thessoloniki,  Greece 

H,N/S 

3214 

17 

1 

2.37 

3 

Honolulu,  Hawaii,  to  Kodiak,  Alaska 

M,N/S 

4057 

85 

3 

2.44 

4 

Davis,  CA,  to  La  Posia,  CA 

M,N/S 

802 

67 

3 

2.46 

5 

Toulouse,  France,  to  Keflavik,  Iceland* 

H,  N/S 

2806 

17 

7 

2.47 

6 

Honolulu,  Hawaii  to  La  Posta,  CA* 

M,E/W 

4287 

20 

7 

2.70 

7 

Coco  Solo,  Canal  Zone  to  Stockbridge, 

NY 

LO,  N/S 

3764 

65 

4 

2.80 

8 

Guam  to  Honolulu,  Hawaii* 

LO,  E/W 

6100 

85,20 

12 

3.01 

9 

Andoya,  Norway,  to  New  Delhi,  India 

M,E/W 

5958 

12 

3 

3.12 

10 

Palo  Alto,  CA,  to  Thule,  Greenland 

TA,  N/S 

5068 

17 

1 

3.18 

11 

Honolulu,  Hawaii,  to  Washington  DC 

M,E/W 

7808 

17 

2 

3.41 

12 

Tarlac,  Philippines,  to  Yokohama,  Japan 

LO,  E/W 

2937 

107 

1 

3.56 

13 

Davis,  CA,  to  Honolulu,  Hawaii 

M,E/W 

3938 

17 

") 

3.58 

14 

Honolulu,  Hawaii,  to  Corona,  CA 

M,E/W 

4200 

95 

1 1 

3.63 

15 

Guam  to  Kodiak,  Alaska 

M,N/S 

7139 

85 

3 

3.80 

16 

Toulouse,  France,  to  Neimakri,  Greece 

M,E/W 

1922 

20 

2 

4.05 

17 

Ft  Monmouth,  NJ,  to  Palo  Alto,  CA 

M,E/W 

4130 

35 

3 

4.06 

18 

Boulder,  CO,  to  Pt  Barrow,  Alaska 

H,  N/S 

481 

125 

-> 

4.29 

19 

Honolulu,  Hawaii,  to  Yokohama,  Japan 

M,E/W 

6193 

1 60, 20 

4 

4.32 

20 

Tarlac,  Philippines,  to  HEH,  Austrailia 

TE,  N/S 

4244 

100 

8 

4.45 

21 

Guam  to  HEH,  Australia 

TE 

5196 

85 

9 

4.50 

22 

Davis,  CA,  to  Kodiak.  Alaska 

M,N/S 

3112 

100 

5 

5.17 

23 

Guam  to  Yokohama,  Japan 

M,N/S 

2505 

85,20 

6 

5.82 

Type  of  path  based  on  location  of  control  points 

M = Midlatitude  All  cases  rms  error  = 3.83  MHz 

E/W  = East,  west 

N/S  = North,  south  * The  data  originally  used  empirically  to  derive  the 

TE  = Transequatorial  MINIMUF  model. 

LO  = Low  latitude  **  1 km  = 0.621371 2 statute  mile  or  0.5399.S68 

H = High  latitude  nautical  mile. 

TA  = Transauroral 
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I iilK  1)1  llic  tO!its  produced  iin  error  r)!'  less  than  3.0  MU/,,  wliicli  is  quite 

eneouraitini!  considering  tlie  variability  ol'day-to-day  MUI's  and  the  makeup  of  the  data 


in  this  set  ( paths  I X ). 

f igure  4 shows  distribution  of  mis  error  for  the  146  path  months  of  data.  This 
histogram  shows  that  the  overall  accuracy  characteristics  of  MlNlMUf'-3  are  centered 
within  the  desired  range  of  1.5  4.0  MH/..  Considering  all  the  variables  intlueneing  the 
■MUf,  the  distribution  eurve  is  quite  well  behaved. 


ALL  CASES 
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CLASS 

RMS 

FREQ 

% FREQ 

1 

0 5-0.99 

0.5 

2 

1-1.49 

7 

4.0 

3 

1.5-1.99 

16 

8.0 

4 

2-2.49 

31 

16.0 

5 

2.5-2.99 

22 

11.0 

6 

3-3.49 

29 

15.0 

7 

3.5-3.99 

32 

16.0 

8 

4-4.49 

19 

10.0 

9 

4.5-4.99 

1 

6.0 

10 

5-5.49 

10 

5.0 

n 

5.5-5.99 

6 

3.0 

12 

6-6,49 

5 

2.5 

13 

6.5-6.99 

0.5 

14 

7-7.49 

4 

2.0 

15 

7.5-7.99 

1 

0.5 

16 

8-8.49 

0 

0.0 

17 

8.5-8.99 

0 

0.0 

18 

9-9.49 

0 

0.0 

19 

9.5-9.99 

0.0 

20 

10-10.49 

0 

0.0 

21 

10.5-10.99 

0 

0.0 

22 

11.0-11.49 

1 

0,5 

RMS  vs  PATH  MONTHS 


IZ2 


1 2 3 4 5 6 7 8 9 10  11  12  13  14  15  16  17  18  19  20  21  22 


RMS  (CLASS) 


Figure  4.  Distribution  of  rnis  error  by  number  of  path  months. 


DATA  TYIM 


A critical  part  of  any  investigation  involving  the  use  of  observed  measurements  is 
the  r|uality  and  time  resolution  of  the  measurements.  This  is  particularly  important  when 
multiple  samples  are  merged  into  mean  values,  as  was  the  ease  with  the  oblique  sounder 
ilata.  As  discussed  in  the  section  on  data  preparation,  there  were  three  types  of  sounder 
data  userl:  ( I ) N I ,SS  I II  DR  ; ( 2 ) N I'SS-stripchart ; and  (3)  non-N  TSS.  To  determine  whether 
data  quality  had  any  impact  on  the  overall  accuracy  results,  table  4 was  prepared  showing 
rnis  error  as  a function  of  sounder  measurement  method. 
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Taole  4.  MlNIMUF''-3  niis  error  versus  sounder  data  type. 


Measurement  Method 

Sample,  path  months 

% of  Sample 

rms  Error,  MHz 

NTSS-HFDR 

59 

30 

3.31 

NTSS-Stripchart 

101 

52 

4.37 

Non-NTSS 

36 

18 

3.27 

As  would  be  expected,  MlNlMUl'-3  provides  better  results  when  compared  to 
better  quality  data.  The  digital  NTSS-HFDR  System  is  automatic  and  can  be  processed 
digitally.  In  general,  the  sounder  instrumentation  is  of  good  quality,  and  care  was  ex- 
ercised to  maintain  the  accuracy  of  the  system.  A path  hour  mean  MUF  for  one  month 
from  this  data  is  a result  of  1800  observations.  While  the  sounding  instrumentation  of 
the  NTSS-stripchart  System  was  of  good  quality,  the  collection  and  reduction  ot  data 
were  loosely  controlled.  In  general,  this  type  of  data  collection  was  an  adjunct  to  normal 
daily  operations  at  a Navy  Communication  Station  and  was  not  monitored  with  respect 
to  a specific  investigation. 

The  non-NTSS  data  were  usually  project-related.  While  the  instrumentation  was  of 
lesser  quality  than  that  employed  with  the  NTSS,  the  close  control  of  the  project  investi- 
gator generally  assured  its  validity.  In  general,  table  4 mflects  this  differentiation. 

3.  PATH  LENGTH 

Figure  5 shows  the  distribution  of  MlNlMUF-3  nns  error  as  a function  of  path 
length.  Between  ranges  of  1000  and  8000  km,  the  rms  error  varies  around  the  desired 
4.0-MHz  accuracy  goal,  ranging  from  an  error  of  4.41  MHz  at  2000  km  to  3.38  MHz 
at  6000  km.  While  fig  5 indicates  that  there  are  no  drastic  range  dependencies  in 
MINlMUF-3,  the  model  appears  to  improve  in  accuracy  with  range.  This  is  probably 
because  there  are  fewer  likely  modes  of  propagation  at  long  ranges  and  the  MINIMUF 
mode  selection  logic  is  more  apt  to  be  right. 

4.  PATH  ORIENTATION 

Table  5 summarizes  the  performance  of  MINIMUF-3  as  a function  of  orientation. 
This  categorization  is  important  to  assure  that  the  sunrise/sunset  reactions  are  correct 
for  varying  degrees  of  path  illumination.  The  north-south  (N-S)  paths  are  those  which 
lie  within  ±15°  of  a 0°  or  180°  bearing.  The  east-west  (E-W)  paths  are  those  which  fall 
within  ±15°  of  a 90°  or  270°  bearing.  The  paths  which  did  not  meet  either  criterion 
were  put  in  the  “other”  category. 

fable  5 indicates  that  MlNIMUF-3  is  slightly  more  accurate  in  predicting  E-W  patlis 
than  N-S  patlis.  fhis  is  not  surprising,  considering  the  dynamics  of  the  abrupt  E-region 
changes  which  occur  when  tlie  entire  path  is  illuminated  suddenly,  as  on  N-S  paths.  How- 
ever, the  magnitude  of  the  error  is  small  enough  to  assume  that  MINIMUF-3  will  produce 
consistent  results  irrespective  of  path  orientation. 
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ALL  CASES  •• 


8 


length 

FREQ 

% FREQ 

RMS 

L < 1 000 

6 

3 

2.46 

1000  <L  <2000 

4 

2 

4.05 

2000  < L < 3000 

28 

14 

4.41 

3000  < L < 4000 

23 

12 

4.08 

4000  <L  <5000 

67 

34 

3.66 

5000  < L < 6000 

26 

13 

4.13 

6000  < L < 7000 

32 

16 

3.38 

7000  <L  <8000 

10 

5 

3.65 

PATH  LENGTH 

••NOTE:  PATH  LENGTH  IN  THOUSAND  km 

Figures.  WINIMUF-.3  Rms  error  as  a function  of  path  length. 


Table  5.  MINIMUF-3  rms  error  versus  path  orientation. 


! Path  Orientation 

Sample,  Path  Months 

% of  Sample 

mis  FTror,  MFlz 

1 Norlh/south 

46 

24 

4.37 

1 Fasi/west 

44 

22 

3.34 

! Other 

1 

106 

54 

3.76 

5.  si:as(3n;month 

F'igure  6 shows  tlie  distribution  of  rms  error  as  a function  of  season  and  month. 
\IINIMUI  -3  is  seen  to  be  slightly  more  accurate  in  the  winter  and  summer  than  during 
the  spring  and  fall.  The  monthly  histogram  in  fig  6 shows  a double  error  peak  at  or  just 
alter  the  eipiinoxes.  The  month  of  October  was  consistently  the  most  difficult  to  predict 
lor  throughout  the  data  set,  which  was  probably  a refiection  of  the  complexity  of  the 
underlying  physical  processes. 
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Figure  6.  MINIMUF  - 3 Rins  error  as  a function  of  season  and  month. 


6.  latitudl 

Tlic  next  tests  were  made  to  determine  error  as  a function  of  geomagnetic  latitude. 
Tlie  three  categories  denote  transequatorial  propagation,  midlatitude  propagation,  and 
transpolar  propagation.  These  general  areas  have  entirely  different  propagation  charac- 
teristics and  problems.  Of  the  three,  the  midlatitude  case  should  be  the  best  behaved  and 
the  most  predictable.  Table  6 indicates  that  this  is  the  case. 


Table  6.  MlNlMUF'-3  nns  error  versus  geomagnetic  latitude. 


(icomagnetic  Latitude 

Sample.  F’ath  Months 

% of  Sample 

mis  Frror,  Mil/. 

25®S  to  2.S*N 
(transequatorial) 

87 

44 

4.32 

20°N  to  60°N 
(midlatitude) 

107 

55 

3.44 

Cireater  than  60°N 
(polar) 

■) 

1 

3.18 
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For  polar  regions,  llie  data  available  were  insuffieient  to  be  eoneliisive.  Furthennore, 
there  were  no  data  from  the  southern  polar  regions,  whieh  probably  would  have  been  quite 
dillerent  than  the  north  polar  region  propagation.  The  most  surprising  aspeet  of  table  6 
is  the  sueeess  MINIMUF'  had  with  transeciuatorial  paths.  The  variability  of  the  F-region  and 
the  eontinuous  l{-region  intervention  around  the  equator  make  those  regions  difficult  to 
prediet.  As  seen  in  table  6,  MlNIMUF-3  does  quite  well  with  these  types  of  paths,  and  the 
data  sample  is  sufficient  to  be  eonelusive. 

7.  SOLAR  SUNSPOT  NUMBFR  (SSN) 

A major  consideration  in  MUF  prediction  is  the  ability  of  the  model  to  deal  with 
dillerent  phases  of  the  solar  sunspot  cycle.  Ideally,  it  should  produce  consistent  results 
lor  SSN  values  between  1 and  150.  In  the  acciuisition  of  data  for  the  MINlMUF-3  verifi- 
cation, care  was  exercised  to  assure  that  a sufficient  cross  section  of  sunspot  numbers 
existed.  Table  7 summarizes  MINIMUF-3  performance  as  a function  of  SSN. 

Table  7.  MlNIMUF-3  rms  error  versus  solar  sunspot  number. 


Sunspot  Number  (Cycle  Phase) 

Sample,  Path  Months 

% of  Sample 

rms  Error,  MHz 

10-30 

(minimum ) 

73 

37 

3.20 

31-60 

(rise  and  decline) 

1 1 

6 

4.15 

61-00 

(near  maximum) 

23 

12 

3.34 

91-120 

(maximum) 

87 

44 

4.34 

121-150 

(high  maximum ) 

2 

1 

4.63 

1 

I able  7 indicates  that  MlNlMUF'-3  accuracy  declines  somewhat  at  the  solar  maximum 

period. 

The  two  path  months  of  data  taken  near  high  maximum  (ie,  1960;  cycle  19)  are 
suspect  because  available  sounder  instrumentation  was  extremely  limited  in  that  period. 

The  remaining  results  are  consistent  with  theory.  As  the  solar  cycle  approaches  max- 
imum, the  level  of  day-to-day  activity  increases.  Compared  to  solar  minimum  conditions, 
solar  maximum  day-to-day  variability  in  emission  and  resultant  ionospheric  support  at 
higher  frequencies  is  larger  and  more  difficult  to  represent  by  mean  values.  This  can  be 
best  illustrated  by  looking  at  the  distribution  of  daily  MOFs  used  to  generate  a mean  MOF 
for  a given  hour.  Figure  7 shows  the  MOF  histogram  for  the  Honolulu,  Hawaii,  to 
( orona.  California,  path  at  1800  UT  for  the  month  of  December  1968.  The  mean  MOF 
is  27.3  MHz  and  the  mis  error  in  the  observed  MOF  is  3.1 1 MHz.  From  the  manner  in 
whieh  the  daily  MOFs  are  typically  distributed,  as  shown  in  fig  7,  it  is  expected  that 
simple  empirical  models  such  as  MINlMUF-3  will  become  less  accurate  at  solar  maximum 
because  ambient  conditions  are  more  difficult  to  typify. 
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Figure  7.  Distribution  of  daily  MOFs  (Honolulu,  HI,  to  Corona,  CA, 
December  1968  - 1800  UT  - SSN  = 98). 


8.  DIURNAL  TRHNDS 

The  last,  and  probably  the  most  important,  test  of  the  MlNIMUF-3  model  involved 
the  study  of  rms  error  as  a function  of  time.  The  objective  was  to  assure  that  the  model 
could  track  through  all  aspects  of  the  day/night  cycle.  To  do  this,  the  entire  data  set  was 
converted  to  local  path  time  (LPT);  ie,  the  local  time  of  the  path  midpoint.  This  put  all  the 
day,  night,  and  transition  characteristics  in  the  same  frame  of  reference.  Figure  8 shows 
a plot  of  rms  error  as  a function  of  LPT.  The  rms  error  varies  between  3.25  MHz  (0400  LPT) 
and  4.6  MHz  ( 1800  LPT).  The  daylight  error  remains  a nearly  constant  3.9  MHz,  with  the 
greatest  variation  occurrring  at  night. 


Figure  8.  MlNIMlJF-3  Rms  error  as  a function  of  local  path  time. 
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I rtic  source  of  the  error  seen  in  the  period  between  1800  and  2400  LPT  is  the  post 

I sunset  lag  in  the  MUF  decay.  This  phenomenon  is  seen  on  most  east-west  midlatitude  paths 

I and  is  a dominant  feature  of  transequatorial  paths. 

1 It  is  most  dittlcult  to  model  such  characteristics  because  they  are  driven  by  ionospher- 

ic transport  processes  which  are  too  complex  for  a simplified  model  such  as  MINlMUF-3  ' 
However,  the  model  does  do  a good  job  of  tracking  the  systematic  changes  which  vary  the 
day, night  MUl-. 

COMPARISON  BLTWFFN  MINIMUF-3  AND  LARGE-SCALE  hf  PREDICTION  CODES 

f 

I Once  the  baseline  accuracy  of  MlNIMLIF'-3  was  detennined  to  be  3.83  MHz  rms 

' error,  the  next  t|ucstion  to  be  resolved  was,  “How  does  this  compare  to  other  hf  MUF 

prediction  tcclmiiiues?”  This  was  a relatively  straightforward  task  to  undertake  because 
I of  a previous  study  performed  at  NOSC.  In  this  study  (ref  17)  three  hf  prediction  large- 

scale  computer  codes  (LSCC)  were  compared  against  obliciue  sounder  data.  Because 
many  of  the  sounder  data  used  were  also  part  of  the  MlNlMUF-3  veritication  data  base, 

1 side-by-side  comparisons  could  be  made.  All  that  was  required  was  to  rerun  MINlMUF-3 

I against  only  that  sounder  infomiation  used  in  the  earlier  study. 

! rile  tlirec  LSCC  programs  reviewed  were  F.SSA-ITSA-I  (ref  3),  which  is  the  basis 

tor  tile  Navy’s  NTP-6.  and  two  versions  of  lTS-78  (ref  2)  designated  the  “red  deck”  and 
“blue  deck”.  These  programs  are  quite  extensive  and  predict  a variety  of  other  propa- 
gation parameters  needed  by  system  designers  and  engineers.  Because  of  their  complexity, 
they  re(|uire  a full-size  computer  (eg,  IBM-370,  UNIVAC  1110,  etc)  to  perform  the  cal- 
culations. 

fo  provide  the  side-by-side  comparison  between  MlNlMUF-3  and  the  three  large- 
scale  programs,  a common  oblique  sounder  data  base  was  identified.  Then  MlNlMUF-3 
and  the  three  large-scale  programs  were  analyzed  in  the  same  manner  presented  in  section 
B.  I he  result  is  compiled  in  table  8 and  is  based  on  137  path  months  of  oblique  sounder 
tlata.  The  data  contain  none  of  the  high-quality  NTSS-HFDR  data,  as  the  earlier  study 
! was  performed  before  this  digitized  sounder  information  was  available. 

Overall,  table  8 indicates  that  MIN1MUF--3  is  comparable  to  the  larger  programs  in 
predicting  MUF.  For  all  cases  tested,  there  is  really  very  little  to  choose  from  among  the 
four  programs.  In  addition,  the  comparison  shown  in  table  8 also  reinforces  a point  made 
in  the  earlier  study  (ref  17 ).  This  was  that,  for  any  given  propagation  situation,  one 
computer  program  will  be  superior  to  the  others.  However,  for  all  situations,  no  program 
will  be  superior  in  all  categories. 

From  table  8,  it  aiqiears  that  MlNIMUF'-3  performs  best  during  uncomplicated 
|)ro|iagation  conditions.  When  geophysical  conditions  arise  which  add  complexity  to  the 
propagation  picture,  the  larger  programs  are  generally  better.  This  is  as  it  should  be 
because  MlNlMUF-3  is  a simplified  estimation  and  is  not  designed  to  handle  every  con- 
ceivable situation.  The  most  important  point  is  that  throughout  table  8,  the  differences 
among  the  four  prediction  programs  in  each  category  are  really  insignificant  with  respect 
I to  solving  day-to-tlay  hf  problems.  File  most  obvious  conclusion  from  this  test  is  that, 

! tor  any  given  protvagation  situation,  a MINlMUF'-3  MUF  prediction  is  as  likely  to  be 

correct  as  one  of  the  large-scale  i^rcdiction  programs  used  in  this  analysis. 


17.  S.illors,  l)b,  HI  Propagation  Predictions:  Program  Proliferation  in  the  Real  World,  paper 
presented  at  N.SA.  1.^  November  1973. 
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Table  8.  Comparison  of  Large-Scale  hf  Prediction  Codes  and  MlNlMUF-3 


RMS  Residual  Error  (MHz) 

Case  Tested 

Number  of 
Path  Months 

1 

MINlMUF-3 

ESSA- 

ITSA-1 

lTS-78 

“Red” 

lTS-78 

“Blue” 

All  cases 

137 

4.03 

5.2 

4.6 

4.5 

Path  length  < 4000  km 

31 

4.41 

5.7 

6.3 

6.5 

Path  length  > 4000  km 

106 

3.88 

4.1 

4.1 

3.6 

4000  < path  length  < .“iOOO  km 

53 

3.88 

5.4 

4.2 

3.7 

Path  length  > 5000  km 

47 

3.88 

4.0 

3.9 

3.5 

25^S  < geographic  latitude  < 25^’N 

59 

4.59 

3.8 

3.4 

3.8 

20^’N  < geographic  latitude  <60^N 

76 

3.65 

5.4 

5.2 

4.9 

Geographic  latitude  > 60‘^N 

1 

3.18 

6.7 

4.8 

4.9 

SSN  < 50 

24 

3.39 

4.8 

4.1 

4.1 

50<SSN  < 100 

58 

4.01 

5.6 

5.3 

5.0 

SSN>  100 

55 

4.30 

4.7 

4.2 

3.9 

Winter  months 

27 

3.64 

4.7 

5.6 

5.1 

Spring  months 

34 

4.08 

5.4 

4.4 

4.4 

Summer  months 

42 

3.92 

5.4 

3.4 

3.5 

Fall  months 

34 

4.40 

4.8 

5.3 

5.7 

DISCUSSION  OF  RESULTS 

Overall,  the  tests  indicate  that  lVllNlMUF-3  produces  a reasonably  consistent 
product  with  a nominal  accuracy  between  3 and  5 MHz  rms  residual  error.  Except  for 
certain  transauroral  cases,  the  model  appears  able  to  handle  a variety  of  propagation 
situations  without  being  overly  sensitive  to  any  particular  parameter. 

The  purpose  of  this  report  is  to  provide  the  reader  with  a new  propagation  pre- 
diction tool  and  with  enough  infonnation  so  that  he  can  confidently  use  it. 

The  reader  should  also  note  that  there  are  geographical  areas  not  covered  in  this 
study  because  of  the  lack  of  readily  available  sounder  data.  Except  for  the  paths  ter- 
minating in  HEH,  Australia,  the  southern  hemisphere  is  almost  totally  excluded.  Trans- 
atlantic and  Asian  land  mass  data  are  also  lacking.  It  is  hoped  that  as  new  data  become 
available,  these  areas  can  be  checked  and  the  results  reported. 


APPLICATIONS 

MINIMUE-3  was  designed  to  provide  a quick  estimate  of  MUF  between  any  two 
geographical  points  for  a variety  of  cases,  using  simple  computational  tools.  It  was  not 
designed  to  replace  the  large-scale  hf  prediction  codes  in  situations  where  the  computer 
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I'acilities  exist.  .MINIMU!-’s  primary  advantage  lies  with  mobile  propagation  forecast  re- 
cinirements  where  the  transmission  patli  is  continuously  changing  and  limited  computer 
I'acilities  exist.  It  is  these  types  of  application  that  will  be  briefly  discussed. 

All  users  of  the  hf  spectrum  arc  concemed  about  whether  a transmission  made  at  a 
given  fretiuency  will  be  heard  at  some  specified  distant  point.  If  the  operating  frequency 
is  in  the  propagation  bandwidth  between  the  LUF  and  MUF,  then  the  transmission  should 
be  heard.  If  not,  some  adjustment  in  fretiuency  must  be  made.  The  primary  problem  is 
rapid  identification  of  the  MUF  and  LUF'  boundaries.  A simplified  method  to  calculate 
l-UF  lias  existed  at  NOSC  for  several  years  and  has  proven  to  be  (luite  satisfactory.  MlNl- 
MUl  -3  provides  an  estimate  of  the  upper  MUF'  boundary  which  allows  a complete  definition 
of  the  propagation  bandwidth  for  any  given  transmission  path. 

To  illustrate  how  the  model  could  be  used,  sample  outputs  from  a mobile  propagation 
forecast  (FROPFllT)  system  being  developed  around  a Tektronix  405  I desk  top  graphics 
calculator  were  generated.  Assume  that  a transmitter  in  Washington,  DC,  wishes  to  know 
his  propagation  coverage  to  London,  Paris,  Madrid,  Rome,  Tel  Aviv,  Capetown,  Omaha, 

San  Diego,  Honolulu,  Bogota,  and  Rio  de  Janeiro  on  16.8  MHz.  Figure  9 (a-d)  shows 
which  of  these  receivers  can  hear  the  transmission  at  0,  6.  12.  and  18  hours  Universal 
lime  (U7  ).  Figure  9a  indicates  that  at  0/hour/UT,  only  those  paths  to  San  Diego, 

Honolulu,  Bogota,  and  Rio  are  open  on  16.8  MHz.  The  remaining  stations  are  in  the  night 
sector  and  the  MUF-  on  each  path  has  dropped  below  the  operating  frequency.  At  06/hours/ 
UT  (fig  9b),  none  of  the  paths  is  open,  since  all  are  in  darkness  or  presunrise.  At  12/hours/ 
UT,  fig  9e  shows  a variety  of  propagation  bandwidth  configurations.  All  the  stations  to 
the  east  of  Washington  are  in  daylight,  and  the  MUF'  is  well  above  16.8  MHz.  On  the  long 
path  to  Ca|)etown.  the  LUF  has  risen  to  above  16.8  MHz.  Stations  to  the  west  of  Washing- 
ton are  still  in  night  darkness  and  the  MUF  has  not  started  to  rise.  The  South  American 
receivers  to  the  south  are  in  a sunrise/early  morning  state.  At  1800/hours/UT  (fig  9d), 
all  the  |)aths  except  Omaha  are  in  daylight  and  1 6.8  MHz  is  a good  frec|uency  for  broad 
geographical  coverage. 

While  only  four  samples  are  shown  in  fig  9,  the  mobile  PROPHET  terminal  just  as 
easily  could  have  produced  24  or  48  examples.  If  driven  by  a date/time  clock,  a display 
similar  to  fig  9 could  be  updated  automatically  at  any  desired  interval.  Manipulation  of 
the  transmitter/recciver  configuration  is  also  quite  easy. 

Simplified  propagation  forecasting,  such  as  that  shown  in  fig  9,  has  many  implications 
in  the  use  of  the  hf  spectrum.  The  ability  to  put  this  capability  in  mobile  platfonns  per- 
mits a continuous  knowledge  of  where  a transmission  can  or  cannot  propagate  to,  allows 
better  use  of  available  freciuencies,  and  provides  a basis  to  exploit  propagation  character- 
istics to  advantage  in  tactical  or  covert  communications.  Because  of  the  small  size  and 
simplicity  of  these  hf  forecast  tools,  it  is  also  realistic  to  envision  that  a point-to-point 
hf  prediction  capability  could  be  an  integral  part  of  future  microprocessor-controlled 
tranmitters  and/or  receivers. 


CONCLUSIONS 

MINIMUI  -3  is  designed  to  complement  existing  large-scale  hf  propagation  codes 
when  computation  resources  are  limited  and  large-scale  codes  are  not  feasible  to  execute. 
In  such  circumstances,  MINlMUI  -3  can  be  used  to  generate  first-order  estimates  of  MUF. 
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Figure  0.  Application  examples  of  MlNlMUF-3, 
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Its  intrinsic  simplicity  (80  BASIC  or  90  F'ORTRAN  statements)  notwithstanding, 
the  accuracy  of  MlNlMUF-3,  as  judged  by  rms  residuals  when  compared  to  hf  oblique 
sounder  data,  is  comparable  to  that  of  existing  large-scale  codes,  l-or  example,  using  a 
common  data  base  consisting  of  137  path  months  of  nondigitized  sounder  data,  MlNlMUF-3 
had  an  rms  residual  error  of  4.0  MH:/.,  wliich  can  be  compared  to  the  best  (4.5  MHz  for  lTS-78 
‘’blue  deck”)  and  the  worst  (5.2  MHz  for  ESSA-lTSA-1 ) errors  experienced  with  the  large- 
scale  codes.  When  MlNlMUF'-3  was  tested  against  59  path  months  of  liigh-quality  digitized 
(NTSS  HFDR)  data,  the  RMS  residual  error  was  3.3 1 MHz.  Overall,  from  a total  data 
base  of  196  path  months  of  data  from  23  different  hf  sounder  paths,  residual  errors  were 
distributed  as  shown  in  table  9. 


Table  9.  Distribution  of  MlNlMUF-3  Errors 


MINFMUF'  rms  Residual  Error 

% of  Total  .Sample 

Less  than  3.0  MHz 

38 

3.1  to  4.0  MHz 

23 

4.1  to  5.0  MHz 

28 

5.1  to  6.0  MHz 

1 1 

For  all  cases,  the  residual  error  was  3.83  MHz. 

For  all  paths,  MlNTMUF-3  showed  no  significant  bias  errors  as  a function  of  the 
local  time  of  the  path  control  points,  which  indicates  that  diurnal  effects  are  being  properly 
modeled  in  an  average  sense. 

When  compared  by  means  of  a common  data  base,  MlNlMUF-3  outperformed  all 
large-scale  hf  codes  on  transmission  paths  of  less  than  4000  km.  Although  this  may  be  an 
artifact  of  the  data  base,  it  does  indicate  that  f^F-i  is  probably  being  modeled  better  than 
the  M factor.  This  is  what  one  would  expect,  considering  the  overly  simple  M factor  model 
employed  in  MINlMUF-3.  (llalso  suggests  that  refinements  of  MlNlMUF-3  should  begin 
with  refinements  in  the  M-factor  algorithms.) 

From  the  overall  data  base,  MlNlMUF-3  accuracy  shows  no  significant  dependence 
on  latitude,  orientation,  or  seasonal  factors.  It  will  degrade  slightly  at  high  sunspot  num- 
bers (ie,  SSN  100  expected  in  the  1980-1982  time  frame)  when  day-to-day  ionospheric 
variability  is  high  because  of  solar  activity.  MlNlMUF-3  accounts  only  for  F-region  pro- 
pagation and  its  predictions  will  degrade  in  the  presence  of  mixed  modes  caused  by 
sporadic  E*  or  when  the  E-layer  itself  dominates. 

MlNlMUF-3  is  based  on  the  simple  idea  that  foF2  'tiodeled  to  a first 

approximation  as  the  lagged  response  to  a driving  function  proportional  to  (cos  x)'^ 
where  x the  instantaneous  solar  zenith  angle  and  where  the  daytime  lag  is  quite  season- 
ally dependent.  The  general  success  of  MlNlMUF-3  over  a broad  range  of  hf  transmission 
path  conditions  lends  some  credence  to  tliese  ideas.  However,  consistent  deviations  from 
single  exponential  decays  at  night  indicate  that  more  complex  relaxation  processes  are 
involved  and  will  be  reviewed  as  future  refinements. 

It  is  concluded  that  MINIMUF-3  is  sufficiently  accurate  to  provide  a simplified 
prediction  of  hf  maximum  usable  frequency  (MUF)  suitable  for  use  on  small  mobile 
propagation  forecast  (PROPHET)  tenninals.  With  its  inherent  simplicity,  it  is  fully  expected 

Because  of  the  well  localized  height  of  sporadic  H reflections,  it  is  a relatively  straightforward  problem 
to  extend  MIN1MUF--3  to  generate  MUF  forecasts  for  a variety  of  mixed-mode  structures. 
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that  other  investigators  will  refine  and  improve  the  model  presented  here.  It  is  hoped  that 
the  successes  along  this  line  will  be  openly  reported  so  that  the  hf  community  as  a whole 
can  benefit. 


RECOMMENDATIONS 


It  is  recommended  that  work  continue  to  further  refine  the  MINlMUF-3  model  to 
include  1 and  sporadic  li  region  propagation.  With  the  acceptance  of  simplified  propagation 
models,  on-site,  real-time  propagation  forecasting  on  a system-by-system  basis  becomes 
realistic  and  cost  effective.  It  is  recommended  that  a microprocessor-based  propagation 
forecast  (I’ROPHI-T)  cajiability  be  considered  as  an  integral  part  of  new  hf  system  develop- 
ment. further,  it  is  recommended  that  work  continue  in  developing  applications  for  the 
new  hf  prediction  capabilities  resulting  from  the  availability  of  simplified  models  such  as 
MlNlMUf-3. 
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APPENDIX  A: 


BASIC  PROGRAM  FOR  MINlMLIF-3 
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Input  variables  for  MINlMUF-3  BASIC  program: 


/ ^ 

^ TT  \ 

LI 

- TRANSMITTER  LATITUDB,  radians  1-  - 

<L1<  5) 

W1 

TRANSMITTER  WEST  LONCIITUDE,  radians 

(-  27r<Wl  <2n) 

l: 

RECEIVER  LATITUDE,  radians  1-  - < L2  < 

5) 

VV2 

RECEIVER  WEST  LONGITUDE,  radians  (-  27r  < W2  < 27r) 

M0 

MONTH  (1  <M0<  12) 

1)(> 

DAY  (1  <D6<31) 

T5 

TIME  (UT),  hours  (0.0  <T5  < 24.0) 

OU  I PUT  MUE,  Mil/. 

S8 

SUNSPOT  COEFFICIENT  (250) 

S9 

SUNSPOT  NUMBER 

PI 

3.141593 

PO 

1.570796 

/■ 
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600  REM  : ''INI  ML'P  tl  MOO  3 tTttttttI  ttt*  ««»  »tf  «»!!»»»  M T « 1 »•«  »*» 

605  H3»1.5i> 

610  K’  = ';iNtLl  ■>♦£  IN'  L3  ■'♦COS'  M ^ » CGJ • L2  > tC  Oi > u;  - u 1 
620  C1=mTM' £0R.  1 ’ ■ K T ' -F 0» ■ 1 -£GN' K 7 < '< 

630  K6=H31C1 

635  IF  Kt=  l THEN  t4£< 

636  Kb=l 
640  K5=l  •*  > 

650  J'?=10O 

660  FOP  Fl  = l '2Tt6''  TO  l-l-''2»K6i  STEP  0.  OSi^Si-l  5 

665  IF  K5=l  THEN  670 

666  k5=0.5 

670  IF  W1=>H2  hNP  W1-W2,PI  ThEH  730 

6S0  IF  w:  -ui  hno  w2-i«ir>Pi  then  730 

630  W0=W1 

700  A=P0-L1 

710  C=P0-L2 

720  GO  TO  760 

730  m=P0-L2 

740  C=P0-L1 

750  W0=W2 

760  B=G1 

770  GOSUB  1500 

780  C=D 

730  e=Kl.fGl 

800  GOSUB  1470 

810  C=B 

820  B=D 

830  GOSUB  1500 
840  W0=W0+[< 

850  L0=P0-e 

860  Y1=0.0172*(104(M0-1 >»30.4+D6) 

870  Y2=0.409»COS('Y1> 

880  K8=3.e2»W04l2+0.  1 3» ^ S IN<  YH  + l . 2»SIN( 2»Y1  > '■ 

830  K3»K8-12*<l4SGN<K8-24)  )»SGN<flBS<  K8-24--  :■ 

831  IF  COS(L04Y2^>-0.26  THEN  300 

836  G0=0 

837  M3*SIN(2.5»G1»K5) 

838  M3»1+2,5»M3*SQR<M3> 

833  GO  TO  1215 

300  K3=(-0.264SIM(Y2>tSlH>;L0>>^aOS>  Y2>»COS<L0  +1 .0E-3) 

310  K3=12-ATN<K3-'SQR<ABSa-K3tK3  - V- ♦■24/PI 

320  T=K8-K3 '2-H2*a-SGN>'1.8-F3  -2  ■ ■ ♦SGN<’hBS' '■"  2'  ■ 

330  T4=K&+F-3  -'2-t2*a+SGN<  K8+('3  2-24  > jTSGN' ABS' r . -.4  . 

340  C0=ABSCCOS>  L0+Y2:' > 

380  T3=3.  7TC0T3-.6 
385  IF  T3>0.1  THEH  330 
388  T3=0.  1 

330  M9=SIN(2.5*Gl»K5i 
1000  113=1+2.  StMSTSQP 019  ' 

1030  IF  T4<T  THEN  1060 

1040  IF  (T5-T>»>'T4-T5:'  '0  THEN  1070 

1050  GO  TO  1240 

1060  IF  a5-T4)»<T-T5 ' >0  THEN  1240 

1070  T6=T5+12*<l+SGHa-T5>  >4SGH'  MB3a-T5>  :< 

1090  G9=PI»<T6-T>/K9 
1100  G8»PI»T3/K9 
1110  U=<T-T6>/T9 

1190  G0=C0»<SIH(G3>+G8*(EXPai''-COS'  G3:' > ■/>  l+GS'HTJ  ■ 

1 1 95  G3=C0T < G8» ( EXP ( -K3/T3 >♦  1 ' '« tEXP > (;K3-24  ^ /2  ' ■ 1 +GS«Gt 
1200  IF  G0=>G3  THEH  1215 
1210  G0»G3 

1215  G2=<l+S9/S6'.-«I19»SQPf6+58»SOP<G0  ■ ' 

1220  G2»G2»<l-0. ltEXP(fK3-24)/3<  ■ 

1225  G2=G2»(l  + <l-SGHai)»SGNa2>)»0.  1> 

1 226  G2=G2» < 1 -0 . 1 » d +SGN ' ABS ( S I N*' L0 ) ^-COS ( L0/'  '•  • 

1230  GO  TO  1430 

1240  T6«T5+12»<1+SGN(T4-T5)^»FGN< ABSa4-T5' > 

1250  G8=PI»T9/K3 


12b0  U=>  T4-T6>  '2 
1330  Ul=-K9  T9 

1400  b0=C01;(.G3t>.EXP  ur'  + r'itES'PiU'  <1+G8TG8.;' 

1420  GO  TO  1215 

1430  IF  G2'J9  ^HEN  1450 

1440  J9=G2 

1450  HEl'T  K1 

1460  RETURN 

1470  REM  ttt  COSLOW  *tt 

1490  D = 0CS(,C0S‘'0''«C0S(B)  + SINf  AUSINfBUCOSiO  ■' 

1490  RETURN 

1500  REM  tt4  TANLOW  ttt 
1510  S=>'h  + B + C ■)  2 

1520  IF  C --1.0E-5  AND  ABSi  C -ABS < A-B  - - >1 . 0E-5  THEN  155 
1530  D’0 
1540  RETURN 

1550  IF  M0i  . S-C '•  -1 .0E-5  THEN  1590 
1560  D=FI 
1570  RETURN 

1580  D=2tMTN'  SQR'  9IN>  S-h>»SIN>  S-B  ^ ^ SIN<  S-C  •tSIN'  S'.’ '< 
1590  RETURN 
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APPENDIX  B: 


FORTRAN  PROGRAM  FOR  MINlMUF-3 
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Input  variables  for  MINIMUF'-3  FORTRAN: 


W TRANSMITTFR  LATITUDF,  radians^-  ^ < W < ^ j 

X TRANSMITTFR  WFST  LONGITUDE,  radians  (-27r  < X < 27r) 

Y - RFC'FIVFR  LATITUDF,  radians  ^ < Y < ^ j 
Z RFC'FIVFR  WFST  LONGITUDE,  radians  (-2it  < Z < 27r) 

M MONTH  (1  <M<  12) 

D DAY  ( 1 < D < 3 1 ) 

T5  TIME  (UT),  hours  (0.0<T5  < 24.0) 

J4  CMUF  - OUTPUT  MUF,  MHz 

S8  SUNSPOT  COEFFICIENT  (250) 

SO  SUNSPOT  NUMBER 
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145 

146  r 

147  n 

148 

149  C 

150  C 

151 

152 

153  C 

154 

155 

156 

157  C 

158 

159 

160 
161 
162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177  100 

178 

179 

180 

181  120 
182 

183 

184 

185 

186 

187 

188 

3^80 240. 


SU99'''UT!Mt;  MuriK  w,  X.  Y.  ?,  m,  0.  T5,  J9,  S8,  59  ) 


uaniCAL  SON 


REiL  K7,  <5,  LO.  xS,  k9,  j9,  K9 

REAL  SAV5(24) 


C0MM/N  /RR0CK/  SAVE 

data  pi  /3. 141593/,  PO  /I. 570796/,  PI  /6. 283185/. 

« DED  /57. 29578/ 


K7  = S!^Jf»l)»SI^(Y)+C7S(4)♦C35(Y)«r^S(^-X) 

G1  = ARCH5(  K7  ) 

M3  = 1.59 
AK6  = H3  » G1 

If  ( AK6  .lT . 1.0  ) AK6  = 1.0 

KMTP  = 1 

ir  ( AK6  .GT.  1.0  ) KM'/P  = 2 
K5=l . 0/AK6 

IP  ( K5  .IvF.  1.0  ) K5  = 0.5 
J9  = ino  . 0 

Da  820  K1  = 1.KMC1P 

AK1  = 1 . 0/ f 2 . 0«AK6  ) +FL5!aT  ( Kl-1  ) * ( 0.9999  - 1.0/AK6  ) 

SGA'  = ? .5“=.  X .and.  2 - X .LT.  PI  .t*R. 

# X .GT.  2 .and.  < - 2 .GT.  PI 

IF  ( .N?)T.  SGn  ) Gl?  Tfl  100 
A = PO  - W 

B = PO  - Y 

WO  = X 
G0  T3  120 
C0NTINUE 
A = PO  - Y 

B = PO  - W 

WO  = 2 
Continue 

alpha  = TANLAW(  a,  Gl,  0 ) 

8 = AK1«G1 

LO  = ARC0S(  C0S(  A )»C0S(  5 ) ♦ SIN(  A )*SIN(  P )»C0S(  ALPHA  ) ) 

delta  = TANLAW(  a,  LO,  9 ) 

LO  = PO  - LO 
WO  s WO  ♦ DELTA 

IF  ( wn  ,GT.  PI  ) WO  = WO  - PI 
_X4-iJl-,-aiZ2a.U  0 , 0 La  A I ( AA-44.  « 38 . 4 ♦ 0 ) 


4 

i 
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19C  Y230.4n9«C  1S<Y1  ) 

191  K8  = 3.S?»i'<  0*12. 0 + 0. 13»(S1N{Y1)+1.2*SIN(2.0*Y1)) 

192  290  K3  = K3  - 12.0o(1.0  ♦ SIGNC  K8  - 2-1 . C ) ) *5 1 GN  ( ABSt  K6  - 24.0  ) ) 

193  ir  ( C?3(  LO  ♦ Y2  ) .GT.  -0.26  ) G0  T0  300 

194  GO  = 0.0 

196  M9  : Sl\'(  2.5*G1«K5  ) 

196  M9  = 1.0  ♦ 2.5«M9oSGRT(  *-19  ) 

197  G0  090 

19°  300  K9  = ( -0.26  ♦ SIN(  Y2  )«SIN(  UO  ) )/(  C0S(  Y2  )*C0S(  LO  ) ♦ 

199  n 1.C9-3  ) 

200  K9  = I'^.n  - ATANC  K9/S0RT(  APSC  l.C  - K9«K9  ) ) )»7,6394 

201  310  T:K°--<9/2.r*12.O»(l.O-SIGN(K8-K9/2.O))»SIGN(A0S(K8-K9/2.O)) 

2C9  330  T4r/-,*K9/2  . n-l2  . n»  ( 1 . 0 + S I GV  ( X8  + K9/2 . 0-24 . 0 ) ) * 

20‘«  « c I G.m  ( AGS  ( k9*k9/2 . 0-24 . 0 ) ) 

204  CO  = A9S(  C?S(  LO  ♦ Y2  ) ) 

20  6 360  T9:9.7(»C0*«9,6 

20A  IF  ( T9  .lT.  0.1  ) T9  = 0.1 

2C7  370  ^•9:S1  M(2.6oGl*K5) 

2'’P  '-9  = 1.0  + 2.p»^‘9»SORT("9) 

209  °1C  IF  (T4.L'f.T)  G0TZ  440 


210 

IF  ( 

(T5-T)«(T4-T5) 

.CT  . 

0.0  ) 

Ge 

T0 

450 

211 

G0  T 

610 

212  440 

IF( 

(T6-T<)*(T-T5) 

.GT. 

0.0  ) 

G0 

T0 

610 

21  3 450  T6  ° T6  + 12.0<>(1  ,0  + ?IGMT-T5)  >«SIGN(  ABS{T-T5)  ) 

21  A G9:°  I • ( T6-  T ) /“■9 

21  5 470  GP-^*  I «T9/'r<9 

216  U = ( T - T6  )/T9 

217  U-4miM(u.  120.0) 

21 e U-Av4x1(U.  -120.0) 

219  660  GO=CC»{SIM(G9)+G8*(EYP(U)-C0S(G9) ) )/(1.0*G6*G8) 

220  G3:C^*(G6»(PXP(-K9/T9)+1.0))oF)(P((io9-24.0)/2.0)/(1.0*G8*G8) 

291  IF  ( GO  .LT.  G3  ) GO  = G3 

222  59G  C2  = ( 1.0  ♦ S9/S8  )«M9oSQRT(  6.0  ♦ 58.0<'SQRT{  GO  ) ) 

r-z  : G?.(  1.0  - O.loFXPI  ( K9  - 24.0  )/3.0  ) ) 

224  r?  r g2«(  1.0  ♦ ( 1,0  - SIGN!  W )*S!GN(  Y ) )*Q.l  ) 

225  G2  = G2*(  1.0  - G.l*(  l.C  ♦ SIGN!  ABS(  SIN(  LO  ) ) - C0SI  LO)))) 

226  G0  T3  800 

227  610  T6  = T5*12 . 0»(  1 . 0 + RIC\i(  T4-T6)  )eSIGM(  ABS(T4-T5)  ) 

2’£  G8=P!»T9/K9 

229  U=(T4.T6)/2.0 

230  U = A8IN1(  U.  120.0  ) 

251  U = AmaxK  U.  -120.0  ) 

232  U1=-a9/T9 

233  IJl  : A«'IM(  Ul,  120.0  ) 

234  Ul  = AKAXK  Ul.  -120. C ) 


235 

770 

GO=CC«(G8*(EXP(U1)*1.0) )*eyP(U)/(l.O*G8*G0) 

236 

Gf!  Ti?  590 

237 

800 

CONTINUE 

238 

IF(G2.GT.JV)  GP  TP  635 

239 

810 

J9sG2 

240 

815 

C0^  T I NUE 

241 

LO  = ( K] 

- 1 )*7 

242 

SAVP(LG*1 ) 

= LO 

243 

SAVE  ( Li^*?  ) 

= WO 

244 

SAVE (LC+3 ) 

= Ke 

245 

SAVE(UQ*4 ) 

= K9 

246 

SAVE (LP*5) 

= T 

247 

SAVE{L0*6) 

= T4 

246 

SAVE(L0+7) 

= CO 

249 

820 

C0NTINUE 

250 

SAVE(15)  = 

G1 

251 

SAVE(16)  = 

AK6 

252 

SAVE(17)  = 

K5 

253 

SAVE(ie)  s 

Y3 

254 

SAVE(19)  = 

Y? 

255 

SAVE(20)  = 

D 

256 

SAVE(21)  = 

FL0AT(KHeP) 

257 

SAVE(22)  = 

FL0AT(M) 

256 

SAVE(23)  = 

S9 

259 

SAVE(24)  = 

SIGN(  W )*EIGN( 

Y ) 

260 

SAVE(25)  = 

S1GN(  ABS(  S1N{ 

LO  ) ) - C0S(  LO 

261 

630 

RETURN 

262 

END 

ERR0RS  0 SI?E  C2447  START  C002? 


ri)f  c ■ I 'Vv  T ^ ► i i - ( t 


■ . r ) 


r' 


7('^ 

?PA 

r 

?tb 

c 

?'^6 

r 

2f'7 

r 

?6P 

?fc9 

27r 

271 

272 

1 00 

273 

274 

120 

275 

276 

277 

278 

140 

279 

280 
281 
282 

160 

283 

284 

285 

ERRTIPS 

r:'  A sfhiPicAL  Tf-iCMa?  »'Ith  siofs  a.  -i.  c. 

TA‘j|  A Ar.T|_Kr<;  -rnr  ANTLc  r'PP^^ITE  C P'  PARIANS 

IF  f A * 5 .'•f.  C .Afr. 

» A ♦ c . ' F , - . A‘  n . 

S p * C . F . A ) F :•  T / 1 ^ 0 

WP ! TA ( A , 1 c' ) A,  P.  C 

FFr.  ►'A  T ( /??!-  I '.'V/ L I r AFr,F  IN  TANLA'P,  1P3E15.5/) 

P F T 1)  - 

IF  ( C .GT.  l,rE-5  .ant.  Af^SC  C - ABS(  A - 9 ) ) .GT. 

« G • T S’  1 4 p 

T A A L A 'A  = C . n 
Pf  TU"^  \ 

S = ( A * P ♦ C )/?.(' 

IF  ( APG(  - C ) .GT.  1,0^  ~b  ) GS^  TPl  160 

TAMAw  = 3.1'11593 

RETUFN 

TaalAW  = P.noATAM  SCPTC  S 1 a|  ( G - A )*SI''I(  S - 9 )/ 

« IF1M5-C)»S1N(1)))) 

RETURN) 

Efjn 

0 PIFC  00355  START  00003 


l.OE-5  ) 


48 


?®fc 

rij'' 

CT  1 - 'j  c ! r,.  ( y ! 

7 

C 

r 

ti  » « 

ooaci/PofiotiaiioPOttPPOPttaocoPPcopoaAuocAPaoaPocdaouottftooo 

2Pg 

r 

« 

2cp 

c 

« 

BB.'O R7P  P VE''- MaROT  , ‘^DSUP 

2Q1 

r 

« 

?92 

C 

» 

0’;f,p  '5  r ; 

293 

C 

• 

:<=,TAI\  TUF  ALG'BRAIC  sign  :'f  THE  A('GU'-lf'iT 

0 

29« 

r 

« 

295 

c 

• 

I Ml-' jT : 

296 

c 

• 

Y s Ah.Y  PEAL  Nl'MPEE 

* 

297 

c 

« 

296 

c 

• 

■’l'Td  ]T : 

299 

c 

« 

PI  ON  s -1.0  IP  Y < 0.0 

o 

300 

c 

• 

pIGNa  C.OIPYsO.O 

• 

301 

c 

• 

PI'.N  r *1.0  IP  Y > 0.0 

302 

c 

• 

303 

c 

• 

304 

c 

305 

c 

306 

IP 

( Y .LT.  0.0  ) SIGN  A -1.0 

307 

IP 

(Y.EO.  0,0)  SIGN*  0,0 

308 

IP 

t Y .GT.  0.0  ) SIGN  = *1.0 

309 

RETURN 

310 

ENO 

ERRPBS 

0 9 

PE  rj0''73  start  00001 

INITIAL  DISTRIBUTION 


NAVAL  ELECTRONIC  SYSTEMS  COMMAND 
PME106-3  (DG  ROBBINS) 

NELEX-OOBE  (W  COOK) 

NAVAL  SECURITY  GROUP  COMMAND 
NAVAL  SECURITY  GROUP  HEADQUARTERS 
G-53  (LT  J FITTEN) 

CHIEF  OF  NAVAL  OPERATIONS 
NOP-941F  (HA  FEIGLESON) 

NAVY  ELECTROMAGNETIC  SPECTRUM  CENTER 
NAVAL  COMMUNICATION  UNIT,  WASHINGTON 
WASHINGTON,  DC  20390 
CAPT  JA  MAOIGAN 

INSTITUTE  FOR  TELECOMMUNICATION  SCIENCES 
BOULDER,  CO  80302 

DON  LUCUS  (3) 

CENTER  FOR  NAVAL  ANALYSIS 
1401  WILSON  BLVD 
ARLINGTON,  VA  22209 

CDRLSHOUP  (3) 

MEGATEK  CORPORATION 
1055  SHATTER  STREET 
SAN  DIEGO,  CA  92106 
DR  PH  LEVINE 

MITRE  CORP 

WESTGATE  RESEARCH  PARK 
MC  LEAN,  VA  22101 
DR  ALLAN  SCHNEIDER 

GENERAL  ELECTRIC  CO 
BLDG  4,  ROOM  48,  COURT  STREET  PLANT 
SYRACUSE,  NY  13201 
DR  GEORGE  H MILLMAN 

TRACOR,  INC 
1600  WILSON  BLVD 
ARLINGTON,  VA  22209 
DR  THEODORE  J COHEN 

SOUTHWEST  RESEARCH  INSTITUTE 
8500  CULEBRA  ROAD 
PO  DRAWER  28510 
SAN  ANTONIO,  TX  78284 
WM  SHERRILL 

COMMUNICATIONS  RESEARCH  CENTER 
SHIRLEY  BAY,  BOX  490,  STATION  A 
OTTAWA,  ONTARIO,  K1  N85S,  CANADA 
Ul  CAMPBELL 

JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
JOHNS  HOPKINS  ROAD 
LAUREL,  MD  20810 
TOM  EVANS 

DEFENSE  DOCUMENTATION  CENTER  (12) 


50 


1 


